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Abstract. 

The characteristic function for the joint measurement of the changes of two 
commuting observables upon an external forcing of a quantum system is derived. In 
particular, the statistics of the internal energy, the exchanged heat and the work of 
a quantum system that weakly couples to its environment is determined in terms of 
the energy changes of the system and the environment due to the action of a classical, 
external force on the system. If the system and environment initially are in a canonical 
equilibrium, the work performed on the system is shown to satisfy the Tasaki-Crooks 
theorem and the Jarzynski equality. 
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1. Introduction 

About one decade ago Jarzynski proved a quite remarkable equality that relates free 
energy differences to average exponentiated work done on a thermally insulated system 
that is acted upon by external time dependent forces varying according to a specified 
protocol. This equation, now commonly referred to as the Jarzynski Equality, reads [1]: 

(e-^-) = e-^^^, (1) 

with w the work, AF the free energy difference between a reference equilibrium state 
of the system at the initial temperature (3 with the force values reached at the end of 
the force protocol and the truly thermal initial state. Note that the Jarzynski Equality 
holds irrespectively of whether the system ever reaches this reference equilibrium state. 

Crooks [2] later showed that Eq. (1) results from the following work fluctuation 
theorem: 

Ptf,to{u!) ^ ^-p{AF-w) /2^ 
Pto,tf{-w) 

that relates the probability density function (pdf) of work ptj,fy{w) in the real forward 
process proceeding from the time to until tf to the pdf of negative work 'pto,tj{—w) of 
the mirror image process where the time ideally runs backward. 

An important question that naturally arises is whether Eqs. (1) and (2) keep 
holding in the more realistic situation where the system remains in thermal contact 
with its environment while the forcing protocol is in action. In regard to Eq. (1) a 
positive answer to this question was already given in Ref. [1] on the basis of classical 
arguments. 

Moreover, quantum extensions of Eqs. (1) and (2) were developed, too, first the 
quantum version of the Jarzynsky equality for cyclic processes [3], and shortly after 
the Crooks theorem was demonstrated to hold also for closed quantum systems [4]. A 
specification of the full statistics in terms of the characteristic function, i.e. the Fourier 
transform of the probability density function (pdf) of work was obtained in [5], and a 
proof of the Crooks relation on the basis of this characteristic function was given in Ref. 
[6]. A generalization to arbitrary initial states such as the microcanonical state was 
obtained in Ref. [7]. In the latter case a microcanonical Crooks theorem was proved 
by means of which the change of the thermodynamic entropy can be inferred from 
the statistics of the work [7, 8]. Illustrative examples for differently driven quantum 
harmonic oscillators were discussed in Refs. [9, 10]. An experimental test was proposed 
in Ref. [11]. 

Quantum generalizations of the Jarzynski equality and the Crooks theorem for open 
quantum systems have almost exclusively been studied for systems with Markovian 
dynamics [12, 13, 14, 15]. Therefore, so far, only systems with weak coupling to 
their environments have been considered as weak coupling is implicit in the Markovian 
assumption [16, 17, 18, 19] | . 

I In the strong coupling limit the dynamics of a particle's position is described by a classical 
Smoluchowski equation [20, 21, 22, 23] but this equation does not provide information about the 
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In the present paper we too restrict ourselves to the case of weak interaction between 
the system and its environment. We do so because to our best knowledge presently 
an unambiguous definition of work and heat of a small open quantum system is only 
known in this very case of weak coupling [24, 25]. However, neither the dynamics of 
the system nor of the bath is restricted otherwise. In this way we allow for general 
non-Markovian dynamics imposed by the bath and arbitrary force protocols which are 
neither constrained to be fast nor slow. 

Here we employ the characteristic function approach of Refs. [5, 7] to address 
the question regarding the validity of the Eqs. (1) and (2) in quantum systems 
in weak contact with their environment. Hence the focus of this study is on the 
characteristic function of the joint statistics of simultaneous measurements of system 
and environmental energies, which are commuting observables. Under the assumption 
of weak coupling the changes of these two energies imposed by the action of an 
external force can immediately be related to the changes of internal energy and the 
heat exchanged with the environment, or, equivalently, to work and heat. Our central 
result is the derivation of a fiuctuation theorem of the Tasaki-Crooks type, for the joint 
pdfs of either internal energy and heat or, heat and work [see Eq. (32) below]. As a 
corollary of this theorem, we recover the Jarzynsky equation (1) and the Crooks theorem 
(2) for the marginal pdf of work. This formally proves that the validity of these relations 
keeps holding for quantum systems weakly interacting with an environment. 

Previous investigations of the Jarzynski equality for open Markovian quantum 
systems [13, 15] are based on the limit of infinitely many repeated measurements of the 
environment's energy. In contrast, the present approach relies on energy measurements 
of the system and the environment at the beginning and the end of the protocol. We note 
that other fiuctuation theorems are known in literature for the heat exchange between 
two systems which either couple directly to each other [26], or through an intermediate 
system giving rise to a heat fiow [27]. In contrast, here we consider systems which 
initially are in thermal equilibrium with their environment and are then driven out of 
equilibrium by a classical external force. 

The paper is organized as follows. In Sec. 2 we define and obtain the general 
expression for the joint pdf of internal energy and heat. In Sec. 3 the joint pdf is 
further evaluated for the initial canonical state of the total system and a fiuctuation 
theorem for work and heat is obtained. Conclusions are drawn in Sec. 4. 

2. Characteristic functions of work, heat and internal energy 

We consider a system 5" that is in weak contact with its environment B. Accordingly, 
the Hamiltonian H{t) of the total system consists of a system and an environmental 
part H^{t) and H^, respectively, and the interaction Hamiltonian H^^, i.e. 

H{t) = H-'^it) + H^^ + , (3) 

dynamics of the momentum. 
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where the interaction is small compared to both the system and the environmental 
Hamiltonians. We assume that the gauge is fixed in such a way that the Hamiltonian 
H^{t) coincides with the system energy despite its time dependence [28]. This 
time dependence is caused by an external change of system parameters according 
to a prescribed protocol. Since the system and the environmental Hamiltonians 
commute with each other the energies of the system and of the bath can be measured 
simultaneously with possible results ef{t) and ef which are the eigenvalues of the 
operators H^{t) and , respectively. The corresponding projection operators onto the 
common eigenfunctions of these operators are denoted by Pi a(t). Hence, the eigenvalues 
and eigenprojection operators are determined by the following equations: 

Hs{t)P,^^{t) = ef (t)P,,„(t), HBP^,a{t) = eaP^At) ■ (4) 

The projection operators onto the common eigenspaces of two commuting hermitian 
operators are hermitian: 

PUt) = P^At) , (5) 

idempotent and mutually orthogonal: 

Pi,a{^)Pi',a'it) = Si^iiSa,a'Pi,a{^) , (6) 

and complete: 

J2P^At) = ^. (7) 

i,a 

where 1 denotes the unit operator on the total Hilbert space of system and environment. 
The first measurement is performed at the time to at which the protocol starts to act. 
At this very time the state of the total system is assumed to be given by the density 
matrix p{to) which we will specify later. The joint probability to measure the respective 
system and environmental energies ef (to) and , in this state is given by 

Pi^a = TrPiAto)p{to) , (8) 

where Tr denotes the trace over the total Hilbert space of the system and the 
environment. After a measurement with the outcome ef (to), the system is found in 
the initial state projected unto the corresponding subspace. This state is given by 

Pi,a = Pj>^i,a(to)p(to)Pi,«(to) . (9) 

The second measurement of the system and the environmental energies is performed at 
the end of the protocol at time tf. By then the density matrix of the total system has 
undergone a unitary time evolution in the total Hilbert space to a new state reading 

PiAtf) = Utf,toPi,aUl^t^ . (10) 

The result of the second measurement is characterized by the conditional pdf of finding 
energies ef,{tf) and ef, given that the result of the first measurement were ef{tf) and 
ef . This pdf is given by 

Ptf,^^', «) = TTPi,^a'itf)piA^f) ■ (11) 



Fluctuation theorems in driven open quantum systems 



5 



Consequently, the joint pdf ptf^^Ae^ , Ae^) to measure changes of the system and the 
environmental energies Ae^ and Ae^, respectively, becomes 

Pt„t,{Ae^ Ae^) = ^ 6 {Ae' - (e^itf) - ef (tp))) 

i,i',a,a' 

X 6 {Ac'' - (ef, - ef )) a'|z, . (12) 

At weak coupling between the system and the environment the random quantity Ae*^ 
determines the change of the internal energy of the system. Then the change of 
the environmental energy equals the amount of energy which is exchanged as heat 
Q = —Ae^ with the system. The small contribution of energy that is possibly released 
from or stored in the interaction Hamiltonian is negligibly small for a system weakly 
coupled to its environment. Hence, the joint pdf of internal energy change E = Ae^ 
and heat exchange Q becomes 

Pt„to{E,Q)= M^-(ef(t/)-ef(to)))5(Q + (ef,-ef))pi^,i„(^',a'|2,a)pi,„. (13) 

,a.a' 

The characteristic function Gf^'^^{u,v) provides an equivalent description of the 
statistics of these energy changes. It is given by the Fourier transform of this pdf, 
which, due to the presence of the delta functions, can readily be performed to yield 

Gf;%{u,v) = JdEj dQe''~^^+^Q%,,{E,Q) 

= J2 e<"('='(*^)-'='(*°)))e-<^(^"-^"))p(^',a'K,a)p,,„. (14) 

i,i',a,a' 

This expression can be further simplified by means of the Eqs. (8), (12) and the 
completeness relation (7) into the form of a correlation function, reading 

Gf;%{u, v) = Tre^("^^f(*/)-^l(*/))e-("^^(*o)-^")p(to) (15) 

where the index H denotes the Heisenberg picture of the corresponding operators , i.e. 

H'nitf) = f//,,„i/^(t/)f/t,,o, Hiitf) = Ul,^H^U,^,o ■ (16) 

The density matrix p(to) describes the system immediately after the first measurement. 
It is given by the projection of the initial density matrix p(to) onto the diagonal states 
with respect to the eigenbasis of the measured operators. It hence assumes the form 

p(to) = Yl PiAto)p{to)PiAto) ■ (17) 

In analogy to the case of the measurement of the total energy of an isolated system 
the characteristic function for a joint measurement of two energies is determined by a 
correlation function of the exponential operator (exp [i{uH^(tf) — vH^)\)H{tf) in the 
Heisenberg picture at the time of the second measurement tf and a second exponential 
operator exp \^—i{uH^ {to) — vH^)\ taken at the initial time. The average is performed 
with respect to the initial density matrix projected onto the diagonal elements with 
respect to the joint eigenbasis of the two operators of the first measurement. 
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In passing we note that in an analogous way as for the internal energy and 
the exchanged heat represented by the Hamiltonians of system and environment, 
respectively, the characteristic function of the changes Xi of any mutually commuting 
set of observables can be represented as a correlation function 

where Ui denotes the Fourier variable conjugate to Xj, X\j{tf) the observable X'^{t) in 
the Heisenberg picture and where the density matrix po is given by the projection of 
the initial density matrix po onto the common diagonal basis of the observables Xj(to), 
2 = 1 . . . as in Eq. (17). 

Once the statistical properties of internal energy changes and the exchange of 
heat are known, the marginal distributions of these quantities can equivalently be 
characterized in terms of their respective characteristic functions Gf^. ^^{u) and G^^ ^^{v)] 
i.e., 

G§,o(n) = Gf;£(n,0). 

G^tM = Gf;l{Q.v) . (19) 

Moreover, the internal energy change and the heat exchange determine the work 

w = E-Q (20) 

performed on the system according to the First Law. Correspondingly, the joint 
characteristic function of heat and work, G^^'^^{x, y), is related to that of internal energy 
and heat by 

G?;Zix,y) = Gf;%{y,x-y). (21) 

Then, the marginal characteristic function of work performed on the system becomes 

G-,^^{z) = Gf;%{z,-z). (22) 

So far we have not yet specified the initial density matrix p(to)- In the next section 
we consider the particularly relevant case of a canonical state of the total system at a 
given temperature. 

3. Fluctuation theorem for work and heat 

We now assume that the total system consisting of the considered system and its 
environment is initially, i.e. aX t = to, in a thermodynamical equilibrium at inverse 
temperature P and is consequently described by the Gibbs state 

P/3(to) = ^-^(to)e-^(^'(*°)+^"+^'"), (23) 

where 

Zito) = Tre-'^(^'(*°)+^"+^'") (24) 
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denotes the partition function. In the particular case of weak couphng between the 
system and its environment a perturbation expansion of the density matrix with respect 
to the interaction Hamiltonian yields up to first order 

P/3(to) ^ P%to) - Zs\to)Z^' f rf^'e-(^-^')(^'(*°)+^")<5i/^^e-^'(^'(*°)+^") , (25) 

Jo 

where 

Z^{to) = Tts e-^^''(*o) (26) 
Z^ = Tib e-^^' . (27) 

Here Tts and Tr^ denote the traces over the system and the environmental Hilbert 
spaces, respectively. The operator 6H^^ = H^^ — {H^^)q specifies the deviation of the 
interaction from its expectation value with respect to the factorizing state 

P°(to) = Z,-i(to)^^^e-^(^'(*°)+^") . (28) 

In order to determine the density matrix p{to), Eq. (17), ppito) has to be projected 
onto the diagonal elements with respect to the joint energy eigenbasis of if'^(to) and 
H^. This projection leaves the unperturbed part p^{to) of the initial density matrix 
(25) unchanged. The first order correction on the right hand side of Eq. (25) vanishes 
in all cases when the interaction between system and environment contains only off- 
diagonal terms with respect to the unperturbed energy basis as is the case e.g. for 
the spin-Bose model [29], or the Caldeira-Leggett model [30] §. In all these cases the 
corrections to the factorizing density matrix p*^(to) are at least of second order in to the 
system-environment interaction. Therefore they can safely be neglected in the limit of 
weak coupling such that the diagonal projection of ppito) leads to the factorizing state 
p°(t), i.e. 

p,(to) = p\t,) ^ Z,i(to)^^^e-^(^^(*o)+^") . (29) 

We emphasize that this holds under the conditions of weak interaction between the 
system and its environment. 

With this initial state the characteristic function for the statistics of internal energy 
and heat becomes 

Gf;l{u,v) = Z^i(to)Z^iTr e<"^H(*^)-^^^(*^))e-<"^'(*")-^'^")e-'^(^'(*°)+^") . (30) 

Analogous to the characteristic function of work performed on a closed system initially 
in a canonical state [6], this expression can be continued to an analytic function of 
the variables u and v on stripes Su = {u = u' + iu"\u' G M, < u" < (3} and = 
{v = v' + iv"\v' G M, —j3 < v" < 0} of the complex plane. The proof is based on the 
fact that expiuH^{t) and exp —ivHs are trace class operators which are analytical 

§ If the interaction contains diagonal terms, these can safely be added to the system Hamiltonian, 
which then would be re-defined accordingly. 



Fluctuation theorems in driven open quantum systems 



8 



functions of u and v in the specified strips [31]. Setting u = u + i[3 and v = v + i(3 one 
finds 

z\to)Gf;lM = z'{tf)Gl;%{-u, -v) 

= Z'{tf)Gf^f^X-u + zP, -V - iP) , (31) 

where is the characteristic function of the internal energy change and heat 

transfer for a fictitious process that runs under the action of the time reversed process 
backward in time. The proof is analogous to the case of a closed system initially in a 
canonical state [6]. Applying the inverse Fourier transform with respect to both variable 
u and V one obtains a Tasaki-Crooks type expression [2, 4]. It establishes a connection 
between the joint pdf of internal energy change and exchanged heat for the original 
process to the corresponding quantity of the time reversed process. This fluctuation 
theorem reads accordingly 

Ptf,toiE, Q) ^ Z^{tf) ME-Q) ^ I3(AF-E+Q) (on\ 

where AF = -(3~Hn[Z^{tf)/Z^'^{to)] denotes the free energy difference between the 
system in the canonical state with the parameter values at t/ and in the initial state at 
to. 



3.1. Nonequilibrium free- energy/work relations 

Replacing the internal energy by the work performed on the system we obtain 

Z'it,)Gf;lix,y) = Z■'^itf)G?J^i-x-y,-y + ^/3) , (33) 
or equivalently 



Here p^.'J(,(Q, w) = ptf^to{Q + w, Q) denotes the joint pdf of heat and work corresponding 
to the characteristic function G^'J^(a:, y) defined in Eq. (21). Integrating over all possible 
values of the exchanged heat one obtains the Tasaki-Crooks theorem for an open system 
which couples weakly to its environment. It reads 



w 



. e-^(^^-) , (35) 

where pfj^tS''^) ~ I pf/ZiQy'^) denotes the marginal pdf of work. From the Tasaki- 
Crooks theorem Jarzynski's work theorem follows immediately saying that the mean 
value of the exponentiated work performed on an open system in weak contact with its 
environment coincides with the exponentiated free energy difference [1], i.e., 

{e-f^-^) = e-"^^ . (36) 
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3.2. Nonequilibrium relations for conditional and marginal probabilities 

From Eq. (32) one can obtain nonequilibrium relations for the marginal distribution of 
heat 

pttoiQ) = I dwp?;l{Q,w), (37) 
for the pdf of work w under the condition that the heat Q is measured, i.e.,: 

PtffyXMQ) = Q , , (38) 

and for the pdf of heat Q under the condition that the work w is measured, i.e.: 
A straightforward calculation yields: 

pUi-Q) 



e^^^(e-'^"'|Q) , (40) 



pttM 

where the symbol {-IQ) denotes average over pt^^t^{w\Q) 
On the other hand one finds the following relation 

Ptf,to{QH _ 



(41) 



Pto,tf{-Q\ - w) 

for the conditional pdf of heat Q given the work w, in the forward and backward 
processes. 

3.3. Nonequilibrium equalities for exponentiated internal energy and heat 
By noting that 

(e^'^)Q = / dQpl,,{Qy^ = Glt^HP) (42) 

and 

{e-^^)E = J dEpl,^^{E)e-^^ = Gg,„(2/5), (43) 

and using Eqs. (19) and (30) one finds two nonequilibrium equalities for heat and 
internal energy: 

^ Zs{to)ZB ^ ^ 

In clear contrast to the Jarzynski equality, the average of the exponentiated heat as well 
as that of the exponentiated energy do depend on the details of the protocol, i.e. these 
relations do not depend only on equilibrium properties of the system. 
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4. Conclusions 

We proved that the Tasaki-Crooks fluctuation theorem and the Jarzyinski equality 
hold in open quantum systems that are weakly coupled to their environments. This 
was accomplished by introducing the characteristic function for the probability of joint 
measurements of system and environmental energies, which are commuting observables. 
This characteristic functions was further evaluated for canonical initial states. Our 
approach rests on two basic assumptions: First, the total system made of the system 
of interest and its environment is governed by Hamiltonian dynamics and second, 
the coupling between system and its environment is assumed to be weak. The flrst 
assumption leads to a unitary time evolution of the total system. The second assumption 
has two important consequences: It allows one to determine the internal energy from the 
Hamiltonian of the system and the heat from the Hamiltonian of the environment. The 
work performed on the system then follows by means of the First Law. The contribution 
of the interaction Hamiltonian is neglected in the definitions of these three energies. The 
second consequence of the weak coupling assumption is the factorization of the state 
of the total system immediately after the first measurements of energies into a product 
of the Gibbs states of the system and the environment as if they were uncoupled up 
to second order corrections in the system environment interaction. In order for these 
corrections to be small for weak but nonzero coupling, the temperature of the initial 
state must not be too low. In contrast, we formally allowed for the exact time evolution 
such that fast as well as slow protocols can be described adequately. 

By an inverse Fourier transformation the joint statistics of internal energy and heat 
(defined as the energy ceded to the environment) was shown to obey a Tasaki-Crooks 
type fiuctuation theorem. From this we derived the Tasaki-Crooks fiuctuation theorem 
and the Jarzynski equality, for a quantum system that weakly couples to its environment 
[6]. Further nonequilibrium relations for the marginal and conditional pdfs of work given 
heat and heat given work were derived. 

In the present paper we restricted ourselves to the case where the external forces 
directly infiuence the system Hamiltonian but leave the interaction and the environment 
Hamiltonians unchanged. Our approach can be generalized to more complex situations 
where the forces immediately infiuence the system plus environment. In this general case 
the Tasaki-Crooks theorem, Eq. (34), continues to hold in a slightly modified form: Also 
the change of the free energy of the environment has to be taken into account, i.e. AF 
has to be replaced by AF + AFb, where AFb denotes the change in free energy of the 
environment due to the parameter change of the environmental Hamiltonian. In those 
special cases where only the system and the interaction Hamiltonians are infiuenced by 
the external forcing, the form of the Tasaki-Crooks theorem as given by Eq. (34) stays 
unchanged. As a consequence the Jarzynski equality remains valid for forces changing 
the system and its interaction to the environment but leaves the Hamiltonian of the 
environment unchanged. 

For a system that strongly couples to its environment it frequently is possible to 
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"dress" the system with excitations of its environment such that the dressed system 
only weakly interacts with its properly redefined environment. Quasi-particles in solid 
state physics are typical examples of this type of systems. An external force acting on 
the original system will also be felt by the dressed system. For the work exerted on 
the dressed system and the heat exchanged with the properly redefined environment 
all conditions for the presented theory apply and therefore the Tasaki-Crooks theorem 
and consequently the Jarzynski equality remain valid for such dressed system. Note, 
however, that for the original "bare" system, as for other strongly coupled systems, 
these relations cannot be expected to hold. 
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